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Disappearance of Squeezed Back-to-Back Correlations - a new signal of hadron 
freeze-out from a supercooled Quark Gluon Plasma 

T. Csorgo 1 ' 2 and Sandra S. Padula 1 
1 Institute* de Fisica Teorica - UNESP, 
Rua Pamplona 145, 01405-900 Sao Paulo, SP, Brazil, 
2 MTA KFKI RMKI, H - 1525 Budapest 114, 
P.O.Box 49, Hungary 

We briefly discuss four different possible types of transitions from quark to hadronic matter 
and their characteristic signatures in terms of correlations. We also highlight the effects arising 
from mass modification of hadrons in hot and dense hadronic matter, as well as their quantum 
statistical consequences: the appearance of squeezed quantum states and the associated experimental 
signatures, i.e., the back-to-back correlations of particle - anti-particle pairs. We briefly review the 
theoretical results of these squeezed quanta, generated by in-medium modified masses, starting from 
the first indication of the existence of surprising particle - anti-particle correlations, and ending by 
considering the effects of chiral dynamics on these correlation patterns. A prerequisite for such a 
signature is the experimental verification that these theoretically predicted back-to-back correlation 
of particle anti-particle pairs are, in fact, observable in high energy heavy ion reactions. Therefore, 
the experimental observation of back-to-back correlations in high energy heavy ion reactions would 
be a unique signature, proving the existence of in-medium mass modification of hadronic states. On 
the other hand, their disappearance at some threshold centrality or collision energy would indicate 
that the hadron formation mechanism would have qualitatively changed: asymptotic hadrons above 
such a threshold are not formed from medium modified hadrons anymore, but rather by new degrees 
of freedom characterizing the medium. Furthermore, the disappearance of the squeezed BBC could 
also serve as a signature of a sudden, non-equilibrium hadronization scenario from a supercooled 
quark-gluon plasma phase. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq, 42.50.-p 
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I. INTRODUCTION 



The major goal of the programme of high energy heavy 
ion physics has been to explore the phases of hot and 
dense, strongly interacting matter, more specifically, to 
find a new phase where quarks and gluons are no longer 
confined to their T = bound states, the hadrons. This 
new phase of matter may be referred to as "Quark- 
Gluon Plasma" (QGP) in analogy to the electromagnetic 
plasma, the phase in which electrons and ionized atoms 
appear as new degrees of freedom, as compared to the 
neutral atomic or molecular forms of matter. 

It has been theoretically predicted by a number of au- 
thors 0, 0, d, 0, H, that the emergence of a large num- 
ber of colored degrees of freedom in an expanding QGP 
could be signaled by correlation measurements. However, 
some of the theoretical expectations related to a first or- 
der phase transition, such as a long-lived QGP, are not 
being observed experimentally. 

Nevertheless, we could still ask the question: can a 
phase transition be signaled by correlation measurements 
at all? Fortunately, the answer is positive: In section 
HI! we shall discuss correlation measurements that have 
indeed signaled experimentally the onset of first order 
as well as second order phase transitions, in solid-state 
physics measurements. 



II. PHASE TRANSITIONS AND 
CORRELATIONS 

Can an experimental measurement of the two-particle 
correlation function be utilized to signal a phase tran- 
sition? A clear yes answer has recently been given in 
condensed matter physics by Schellekens, Hoppeler, Per- 
rin, Viana Gomes, Boiron, Aspect and Westbrook 
In their measurement, a cloud of ultra-cold gas of meta- 
stable He atoms was released from its magnetic trap, and 
after 47 cm ballistic fall, the position of each atom was 
detected. Single particle detection of these neutral atoms 
was possible as the 20 eV internal energy of each atom 
was released when the atoms hit the detector. 

The initial temperature has been varied from above to 
below the Bose-Einstein condensation threshold, Tbec — 
0.5 fiK. The observed two-body correlation function for 
initial conditions with T > Tb ec showed a bunching be- 
havior, while the correlation function became flat for the 
coherent, Bose-Einstein condensed sample, with initial 
conditions T < Tbec- The observed quantum statisti- 
cal correlations are the atomic analogue of the Hanbury 
Brown - Twiss effect [l4[ , and in the case of atomic Bose- 
Einstein condensation, the disappearance of the bunch- 
ing behavior signaled the phase transition from the usual 
thermal state to the Bose-Einstein condensate. The in- 
tercept parameter of the correlation function served as 
an order parameter of this phase transition. Thus the 
correlation function can signal a phase transition clearly. 
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III. FOUR POSSIBILITIES FOR TRANSITION 
FROM DECONFINED TO HADRONIC MATTER 



What kind of phase transitions may occur in hot and 
dense, nuclear matter? According to recent lattice QCD 
calculations, the three kind of equilibrium transitions are 
possible on the (T, plane: if the temperature is in- 
creased at zero, or nearly zero chemical potentials, the 
transition from confined to deconfined matter is a cross- 
over (l5j . i.e. this is not a phase transition in the strict, 
thermodynamical sense: quark and gluon degrees of free- 
dom are present below T c , hadrons are present above T c 
and various observables yield different estimates for the 
critical temperature itself. In contrast, at rather large 
bariochemical potentials, the transition from hadronic 
matter to deconfined matter is of first order. The line 
of first order transitions ends at a critical end point 
(CEP) that separates the first order phase transitions 
from the cross-over region. At this critical end point, 
{(J-CEPtTcep) , the transition from confined to decon- 
fined matter is a second order phase transition. There 
are other more exotic states like color superconducting 
quark matter, or there is a nuclear liquid-gas phase tran- 
sition, but the region of the phase diagram that is rele- 
vant for high energy heavy ion collisions at CERN SPS, 
RHIC and LHC, is the region around the CEP 0. Re- 
cent calculations and experimental data indicate that in 
central Au+Au collisions at RHIC, the transition is likely 
a cross-over: The extracted central temperatures for the 
0-5 % most central ^Jsnn — 130 GeV Au+Au data sam- 
ple are T = 214 ± 7 MeV, and (i B = 77 ± 38 MeV, 
respectively [ljj ■ These values are to be contrasted also 
with the location of the critical end point of the line of 
the 1st order QCD phase transition, located preliminar- 
ily at T CE p = 162 ± 2 MeV and hcep = 360 ± 40 MeV, 
according to the lattice QCD estimate (l6j . 

If the evolution of the strongly interacting matter pro- 
duced in relativistic heavy ion collisions happens through 
thermodynamically equilibrated states, first- and second- 
order phase transitions, as well as an analytic cross-over, 
are the only three possibilities for the state of the sys- 
tem. However, heavy ion collisions create violently ex- 
ploding mini-bangs, where the time-scales of the expan- 
sion are rather short as compared to the typically ~ 100 
fm/c nucleation times of hadronic bubbles in a first or- 
der phase transition from QGP to a hadron gas. Hence 
non-equilibrium transitions are also possible. In ref. [ill ], 
another scenario is suggested, where a rapidly expanding 
QGP state might strongly supercool in Au+Au collisions 
at RHIC and then hadronize suddenly while emitting a 
flash of pions. The predicted signals of this scenario are 
not, even currently, in disagreement with RHIC data on 
Au+Au collisions. So the possibility of non-equilibrium 
hadron formation must be kept in mind and further ex- 
perimental tools must be found to pin down if this mech- 
anism is in fact present in Au+Au collisions at RHIC. 



IV. CORRELATION SIGNATURES FOR THE 
FOUR REHADRONIZATION SCENARIOS 

Let us itemize the possible types of rehadronization 
transitions, and summarize their correlation signatures: 

• A strong 1st order QCD phase transition has been 
studied in the predominant fraction of the litera- 
ture. Such a transition is characterized by long 
life-times and time spreads, since the latent heat 
and the initial entropy densities are large. The pic- 
ture behind this scenario is that of a slowly burn- 
ing cylinder of QGP p. In such a scenario, the 
extended duration of particle emission adds up to 
source parameter in the direction of the average 
momentum of the pion pair, resulting in a substan- 
tial increase of the effective source size. This is 
signaled as R out 3> Rside, regardless of the exact 
details of the calculation [11,13, S EL IE @|- How- 
ever, experimentally R ou t ~ Rside both at CERN 
SPS pi, [H and at RHIC 0, H H|. Hence, 
a strong first order transition from deconfined to 
hadronic matter seems to be excluded by current 
correlation measurements in high energy heavy ion 
collisions, both at CERN/SPS and at the RHIC 
energy range. 

However, three other alternatives still remain. 

• A second order QCD phase transition from decon- 
fined to hadronic matter has been considered re- 
cently in ref. [23| . The interesting conclusion is that 
this phase transition is not signaled by the scale pa- 
rameters of the correlation function, as the spatial 
correlations develop a power-law tail in these kind 
of transitions, and power-laws have no characteris- 
tic scale. Instead, the second order phase transi- 
tions are characterized by critical exponents. One 
of these, traditionally denoted by r], characterizes 
the tail of the s pat ial correlations of the order pa- 
rameter. In ref. [23j, this exponent was shown to be 
measurable with the help of the two-particle cor- 
relation function in momentum space, and it was 
shown that rj — a, where a stands for the Levy 
index of stability of the correlation function itself, 
i.e., C(q) = l + Aexp[-(gi?) Q ]. 

A strong decrease of the Levy index of stability 
a down to about 0.50 ± 0.02 in the vicinity of 
the critical end point is a theoretical prediction 
for the localization of this outstanding landmark 
of the QCD phase diagram. This value was ob- 
tained based on universality class arguments. Ra- 
jagopal, Wilczek and others argued that the univer- 
sality class of QCD at the critical point is that of 
the 3-dimensional Ising model [24|. The exponent 
of the correlation function rj, however, is extre mely 
small in that model, i.e., r, = 0.03 ± 0.01 
In a violent heavy ion collisions, random external 
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fields are also present, which change the univer- 
sality class and, thus, might increase the value of 
the correlation exponent. In the 3-dimensional ran- 
dom field Ising model, rj increases to 0.50 ± 0.05. 
When interpreted as a Levy index of stability, it 
still corresponds to an extremely peaked correlation 
function, which should be very clearly measurable. 
Nevertheless, the excitation function of the shape 
parameter a has not yet been determined experi- 
mentally. Therefore, it is possible that at certain 
colliding energies, either at CERN/SPS range or 
at RHIC, there is a dramatic change in the shape 
of the correlation function, which has not yet been 
identified. 

• A cross-over from QGP to confined matter is the 
last remaining possibility of an equilibrium phase 
transition. This scenario can be signaled by emis- 
sion of hadrons from a region above the critical 
temperature, T > T c , or, by finding deconfined 
quarks or gluons at temperatures T < T c . Re- 
cent lattice QCD calculations suggest that, at zero 
baryochemical potential, the transition from quarks 
to hadrons is a rapid cross-over, different observ- 
ables yielding different transition temperatures at 
Hb = 0. The peak of the renormalized chiral sus- 
ceptibility predicts T c =151(3) MeV, whereas criti- 
cal temperatures based on the strange quark num- 
ber susceptibility, and Polyakov loops, result in val- 
ues higher than this by 24(4) MeV and 25(4) MeV, 
respectively. Signs of quarks or gluons below these 
temperatures have not yet been observed exper- 
imentally. However, emission of hadrons from a 
small but very hot region, with T > T c ~ 176 ± 7 
MeV HH, has been suggested by Buda-Lund 
hydro model fits to the identified particle spectra 
and two-pion correlation functions in Au+Au col- 
lisions at RHIC [13 ■ 

• The fourth possibility corresponds to hadron for- 
mation out of thermal equilibrium. A sudden re- 
combination from quarks to hadrons has been con- 
sidered as the mechanism for hadron formation in 
the ALCOR model [26]. Other realizations of sud- 
den hadron formation and freeze-out were tools 
used to explain the observed scaling properties of 
elliptic flow (2?| with the number of constituent 
quarks. In the context of femtoscopy analyses, it 
has been suggested that the comparable magnitude 
of R ou t ~ R-side ~ Riong could be a signature of a 
flash of hadrons (pion-flash) from a deeply super- 
cooled QGP phase If]. If this mechanism would 
indeed be responsible for the hadronization in high 
energy heavy ion collisions, it would also predict 
that the freeze-out distributions of temperature, 
flow, and density would not depend on the parti- 
cle type (or cross-section), since all hadrons would 
be produced in the same flash. Strangeness would 
be enhanced, as strangeness production predom- 



inantly happens in the pre-hadron state of mat- 
ter. Finally, no in-medium mass modification of 
hadrons could be observed, since deeply super- 
cooled QGP is a state of over-stretched matter, has 
negative pressures internally, which cause the sud- 
den break-up and coalescence of matter. However, 
in this picture, the hadron gas is produced in a large 
volume, and re-scattering effects are small. Hence, 
the interactions which would cause the in-medium 
mass modification of hadrons, do not take place, 
as the hadron gas would freeze out as soon as it 
was produced. Note, however, that the thermody- 
namic considerations in Ref. [Tl[ would also al- 
low, in about 50 % of the parameter space, for the 
production of a super-heated hadron gas. There- 
fore, the indication of a hot spot in the central re- 
gion, and of particle emission from a region with 
T{x) > T c [Tjj would be compatible with this non- 
equilibrium hadronization mechanism as well. 

In the end of this review, we shall propose a new 
method to observe the onset of sudden hadronization 
from a supercooled QGP state as a possible rehadroniza- 
tion mechanism. 



V. COHERENT STATES 

In this session, we will briefly present the theory of 
quantum optical coherence and chaos, by introducing 
the concept of coherent states and an explanation of the 
HBT effect. Then we discuss squeezed states and review 
their applications in high energy heavy ion and particle 
physics. At the end we will focus on how these quantum 
statistical correlations of squeezed hadronic states could 
be experimentally used as tools to search for a sudden 
freeze-out of hadrons from a super-cooled QGP state. 

Let us mention an inspiring historical review by Michel 
Martin Nieto [28| , who compared the well known discov- 
ery of the coherent states by E. Schrodingcr in 1926 [2^| 
with the much less well known discovery of squeezed 
states by E. H. Kennard in 1927 (30| . According to 
Nieto [2g|, Schrodinger's original discovery of coherent 
states was inspired by a question from Lorentz, in a let- 
ter on May 27, 1926, in which he lamented the fact that 
Schrodinger's wave functions were stationary, and did not 
display classical motion. On June 6, 1926 Schrodinger 
replied that he had found a system with classical motion, 
and sent Lorentz a draft copy of his paper, published in 
Ref. [29[. According to it, coherent states can be de- 
fined with the help of the so-called minimum-uncertainty 
method. In this, the mean position and mean momen- 
tum are required to follow the trajectory of the classical 
motion corresponding to the same Hamiltonian, and also 
that the vacuum state is a member of the set of states. 

The development of understanding about the dual na- 
ture of light, evident in its wave-like properties and quan- 
tized detections has been summarized recently in ref. [3l| . 
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When the tools to handle quantum electrodynamics 
were developed, they were applied mainly to high en- 
ergy processes, and it was still naively assumed, that the 
conflict between Maxwell's treatment - who focused on 
the wave-like properties of the electromagnetic fields, and 
whose formulas were the basis of radio engineering - and 
Planck's theory which emphasized the quantum nature 
of light, would be of no significance in optical observa- 
tions. This state of blissful indifference [3l[ was how- 
ever changed by the landmark experiment of R. Hanbury 
Brown and R. Q. Twiss [3], who proposed an interfer- 
ometric method to determine the angular extension of 
distant stellar objects. They found out that intensity 
correlation between photocurrents, recorded in two sepa- 
rated detectors, displayed a bump when the difference in 
optical path lengths between the signals was decreased 
towards zero. These observations have shown that the 
photons in the two different detectors were correlated, 
although they stemmed from two different surface ele- 
ments of distant stars. This way, individual photons, as 
well as pairs of such photons, entered the realm of obser- 
vational optics. Ever since, this phenomenon is known as 
the HBT effect. 

The correlated emission of particle pairs is a funda- 
mental property of quantum fields. The effect has been 
observed in particle physics and been explained based on 
the bosonic nature of pions by Goldhaber, Goldhaber, 
Lee, and Pais in 1960 [32| . Recently, such positive cor- 
relations have been observed also in ultracold quantum 
gases by Schellekens et al [lH, as detailed in Section UTT 
which indicates that quantum fluctuations and correla- 
tions are important basic properties of matter fields, jus t 
as well as that of the electromagnetic (photon) field [33| . 

After the discovery of the HBT effect and the discovery 
of lasers, there was a theoretical debate in the literature, 
arguing if photons from a laser light would show the HBT 
effect or not. The correct theory was published first by 
Glauber who had shown that there is no HBT effect in 
the coherent fields of lasers. In fact, the lack of second- 
and higher-order HBT type of intensity correlations de- 
fines quantum optical coherence. Thus, Glauber related 
the bunching properties of the HBT observations to the 
random, chaotic nature of the photon field in thermal 
radiation. He also pointed out that more information is 
necessary to characterize the quantum state of the pho- 
ton field: "Whereas a stationary Gaussian stochastic pro- 
cess is described completely by its frequency-dependent 
power spectrum, a great deal more information in the 
form of amplitude and phase relations between differing 
quantum states my be required to describe a steady light 
beam. Beams of identical spectral distributions may ex- 
hibit altogether different photon correlations or, alterna- 
tively, none at all. There is ultimately no substitute for 
the quantum theory in describing quanta ". Based on 
Glauber's theory, the puzzling HBT effects observed in 
thermal radiation and the lack of HBT effect in lasers 
were understood in the same framework, and the char- 
acterization of quantized electromagnetic radiation has 



been reduced to counting photons. The new field of quan- 
tum optics was born. These days, we witness the birth of 
quantum atomics [33j as coherent and incoherent beams 
of fermionic (34[ or bosonic [l3[ atoms enter the realm 
of experimental investigations. 

In this section, we highlight some of the basic proper- 
ties of coherent states, introduced by Glauber in Ref. (3f| 
as eigenstates of the annihilation operator. For simplic- 
ity, let us consider a given mode of a free boson field, 
corresponding to a harmonic oscillator. After quantiza- 
tion and normal ordering, the Hamiltonian operator of 
the one-mode harmonic oscillator is written as 

H = ua f a, (1) 

where and a are the creation and annihilation opera- 
tors, respectively. Their non-vanishing commutator is 

[a, a t] = l. (2) 

With the help of these creation and annihilation opera- 
tors, a Fock space can be built up by considering that 

«|0> - 0, (3) 

\n) = -±=a!»\0). ( 4 ) 

The creation and annihilation operators step up or down 
on the infinite ladder of Fock spaces as 

a f |n) = V^TT|n+l), (5) 
a\n) = ^Jn \n — I ). (6) 

Here n stands for a non-negative integer number. The 
number operator is the Hermitian N = a^a, with N\n) — 
n\n). These Fock states are orthonormal, 

(n\m) = S n<m , (7) 

where 8 n , m stands for a Kronccker-delta. The resolution 
of the unity operator in terms of Fock states is given by 

oo 

l = X>)(n|. (8) 

So, these Fock states form a complete, orthonormal basis. 
The density matrix of the system can be expanded as 

oo 

P=^2Pn\n)(n\, (9) 
n=0 

where p n stands for the probability that the quantum 
mechanical system is in state |n). Its normalization is 

oo 

Trp = ^p„ = l. (10) 

n=0 

The above formulas can be generalized with the help of 
coherent states, but with certain subtleties. Let us start 
first with definitions. 
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Coherent states are the eigenstates of the annihilation 
operator, 



a\a) = a\a), 



(11) 



where a is a complex (c) number. Their algebraic prop- 
erties are very interesting. For example, one can express 
them in terms of Fock states as 



expHa| 2 /2)£^=|n) 

71 = 



(12) 



It is interesting to observe that, removing (detecting) 
one quantum from a coherent state does not change the 
probability that yet another quantum could be removed 
from such a state. This is in sharp contrast to the prop- 
erties of Fock states where, once a particle is detected, 
the resulting Fock state is orthogonal to the Fock state 
before the detection, i.e., (n\n — 1) = 0. 

The expansion of \a) in terms of Fock states can be 
used to prove that these states are properly normalized, 
(a\a) = 1. Although different coherent states are not 
orthogonal, i.e., 



|H/3)| 2 =exp(-|a-/3| 2 ), 



(13) 



they can, nevertheless, also be used to resolve unity as 

d 2 a 



\a){a\, 



(14) 



where d 2 a = dRe[a]dlm[a]. It then follows that any 
state may be expanded linearly in terms of coherent 
states. Therefore, the most general light beam (of a given 
mode k and of a given polarization) can be described by 
a density operator of the following form: 

p = J d 2 ad 2 a'V 2 (a,a')\a)(a'\. (15) 

In many practical important cases, such as the ther- 
mal or the coherent radiation, the density matrix turned 
out to be in a diagonal representation when expanded in 
terms of the overcomplete set of coherent states: 



d 2 aV(a)\a)(a\ 



(16) 



A similar representation was proposed by Sudarshan [36j 
shortly after Glauber's work, which emphasized the ad- 
vantage of this representation when evaluating expecta- 
tion values for normal ordered products of creation and 
annihilation operators, 

TrpO = Tr{p(a t ) A a^} = [ d 2 aV(a)(a*) x a^ . (17) 



Consequently, the evaluation of expectation values of 
normal ordered operators in a quantum mechanical sys- 
tem, defined with the help of its density matrix, is re- 
duced to the evaluation of expectation values of the quasi- 
probability distribution V(a), defined over the complex 



plane of a. The Hermiticity of the density matrix implies 
that V(a) is a real valued function, although it is not nec- 
essarily positive, due to the quantum nature of the field. 
This representation of the density matrix is frequently 
called as the Glauber-Sudarshan representation. 

Glauber investigated a model of photoionization of a 
pair of atoms, labeled 1 and 2, lying at ri and r 2 within a 
light beam of sufficiently narrow spectral bandwidth and 
given polarization. Summing up the transition proba- 
bilities over final electron energies, there is no quantum 
mechanical restriction in defining the time at which each 
electron emission takes place. The probability density for 
ionization of atom 1 at time t\ and for atom 2 at time t% 
can be written as 



P 2 (ti,t 2 ) =Pi(ii)Pi(i2)C 2 (*i,*2), 



(18) 



where P\(t) is the transition probability of each atom 
placed individually in the beam, and C 2 (ti, t%) is the two- 
particle correlation function. 

Glauber pointed out, that for a coherent state, corre- 
sponding to V{a) = S(a — (3), C 2 (ti,t 2 ) = 1, i-e., there 
is no HBT type of correlation in a coherent beam, which 
would be a reasonable model for the laser field. The den- 
sity matrix of a thermal radiation with mean number of 
photons (n) can be written as 



= ^Pn\n){n\, 



n=0 



Pn 



(l + (n))i+»' 



(19) 
(20) 



This corresponds to a Gaussian distribution in terms of 
the V{a) representation, 



V{a) 



1 



7r(n) 



exp(-|a| 2 /(n». 



(21) 



It is easy to show that for such a filtered, single mode, 
thermal radiation, the correlation function is 



C 2 = 1 + | exp (ico[t 2 -h- (x 2 - xi)/c}) | 2 



(22) 



This corresponds to the classical limit, when modes are 
treated as forming a continuum with a stochastic noise. 
Note that for such a thermal light, the maximum value 
of the correlator is C 2 = 2, in a sharp contrast to the 
C 2 (t\,t 2 ) = 1 value, that is characteristic of a coherent 
light. 



VI. SQUEEZED CORRELATIONS IN PARTICLE 
FEMTOSCOPY 



Squeezed states are generalized coherent states with 
very interesting quantum statistical properties. They 
can be obtained by means of the displacement opera- 
tor method, as discussed below, by using a Bogoliubov- 
Valatin transformation, or even as generalized minimum 
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uncertainty states [28| . Let us consider first the displace- 
ment operator method, and discuss their properties in 
terms of Bogoliubov transformation in the next subsec- 
tion. 



A. SQUEEZED COHERENT STATES 

The displacement operator has the following properties 

D^ifiaDiP) = a + 13, (23) 
D^ 1 (f3)a^ D(f3) = a^+(3*. (24) 

It is easy to show that acting with the first of the above 
equality on a vacuum state, the displacement operator 
D(a) creates a coherent state, i.e., 

D(a)\0) = \a) (25) 

On imposing the requirement that D(a) should be a uni- 
tary operator, this leads to 

D{a) = exp(aa^ — a*a) 

= exp(— |a| 2 /2) exp(aa^) exp(— a*a), (26) 

where the second part expresses the displacement oper- 
ator in a normal order, i.e., where all the creation op- 
erators stand to the left and all the annihilation opera- 
tors stand to the right. Further discussion of the prop- 
erties of these displacement operators, can be found in 
Refs. [13,|33. 

The displacement operator is generalized to the 
squeeze operator S(z) as 

S(z) = exp{-(2aV - z*aa)}, (27) 
z = rexp(i<f>) (28) 

Squeezed states are defined in terms of displacement and 
squeeze operators as follows: 

\a,z) = D(a)S{z)\0), (29) 

where the ordering D(a)S(z) differs from S{z)D{a) by a 
change of parameters, being however equivalent. 



B. Surprizes on n -tv Bose-Einstein correlations 

In 1991, Andrew, Pliimer and Weiner [39l] pointed out 
to the surprising existence of a new quantum statisti- 
cal correlation among boson- antiboson pairs. The sur- 
prise was related to the fact that this type of correlation 
involved particle- antiparticle pairs, differently than the 
better known Bose-Einstein Correlations (BEC), which 
occurred between two identical particles. We can under- 
stand the origin of this effect in a simple way in terms of 
creation and annihilation operators, taking the 7r° case 



for illustration. For instance, the single-inclusive distri- 
bution is written as 

iVi(kj) = w k ,^- = uj ki (al t a ki ) , (30) 

and the two-particle distribution, after the decomposition 
that follows from Wick's theorem, is written as 

N 2 (ki , k 2 ) = Wki Wk 2 (d kl al 2 d k2 d kl ) 
= w kl w kl { (4, Ski ) (4 2 «k 2 ) + (a^ «k 2 ) <4 2 «k! ) 

+ (^k^O^ki 6 ^)} • (31) 

In Eq. (|3~lj) , the first term corresponds to the product 
of the two single- inclusive distributions, the second one 
gives rise to the Bose-Einstein identical particle correla- 
tion, reflecting their last position just before being emit- 
ted. The third term, absent in the ^ir^ case, is the re- 
sponsible for the particle- antiparticle correlation (either 
7r ± 7r T or 7r°7r°, since ir° is its own antiparticle). They 
are related to the expectation value of the annihilation 

(creation) operator, i.e., to < a k ^ki > ^ u j analogous 
to what is observed in two-particle squeezing in optics, 
where the averages are estimated using a density matrix 
that contains squeezed states, as briefly discussed in the 
previous section. Under the conditions as those usually 
considered in femtoscopic analyses, the last term in Eq. 
(|31[) vanishes. However, it is non zero if the Hamiltonian 
of the system is of the type H = Hq + Hi, where Hq is the 
free part in the vacuum, corresponding to final particles, 
and Hi represents the interaction of quasi-particles, re- 
sulting in an effective shift of their masses. Alternatively, 
as in the pioneer work in Ref. [391 ], this could be similar 
to having a chaotic superposition of coherent states and 
a density matrix containing squeezed states. 

Although that initial discussion by Andreev, Pliimer 
and Weiner was not entirely correct, they clearly pointed 
out that the particle-antiparticle correlations would man- 
ifest themselves as an enhancement above unity of the 
correlation function, i.e., C(ir + ir~) > 1 and C(7r°7r°) > 1, 
reflecting particle-antiparticle quantum statistical effects. 

In 1994, Sinvukovfiol]. also discussed a similar effect for 
7r + 7r~ and 7r°7r° pairs, claiming that it would be due to 
inhomogeneities in the system, as opposed to homogene- 
ity regions in HBT, which comes from a hydrodynamical 
description of the system evolution. He used Wick's theo- 
rem for expanding the two-particle inclusive distribution 
in terms of bilinear forms. 

In 1996, Andreev and Weiner [4lj elaborated further 
their original idea. They considered that in high energy 
collisions, a blob of strongly interacting pions, which was 
seen as a liquid, was formed and later suffered a sudden 
breakup, so that the pionic system having ground state 
and pioninc excitations was rapidly converted into free 
pions. In the moment of transition, they postulated that 
the in-medium creation and annihilation operators (b' , b) 
could be related to the corresponding free ones (a^, d), by 
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means of a squeezing transformation 



cosh[r(i,j, x)] and Si.j = sinh[r(i, j,x)], where 



a = 6cosh(r) + w sinh(r) 
a) = Ssinh(r) + W cosh(r) 



(32) 



where r = h\n(Ef r /Ei n ) is a squeezing parameter (42j, 
Ef r and Ei n are the asymptotic (free) energy and the 
in-medium energy, respectively. 

In the same year, Asakawa and Cs6rgo[43| proposed 
a similar structure to this previous approach, but re- 
lating in-medium operators to free ones by means of 
a two-mode Bogoliubov transformation. They also 
proposed to observe hadron mass modification in hot 
medium by means of Back-to-Back Correlations, relat- 
ing particle-antiparticle pair correlations to two-mode 
squeezed states. 

There were a few more tentative works by the two 
groups but the correct approach was finally written 
in 1999, by Asakawa, Csorgo, and Gvulassv 44j . The 
squeezing parameter they proposed was somewhat simi- 
lar to the one described by Eq. (|32|) . 



C. Back-to-back boson-antiboson correlations 

The formalism developed by Asakawa, Csorgo and 
Gyulassy for squeezed bosons in an infinite medium can 
be summarized as follows. The in-medium Hamiltonian, 
H, is written as H = Hq — \ J dx dy(f>(x)5M 2 (x — 

y)0(y), where H = \ J dx (j) 2 + |V0| 2 + m 2 2 ) , is 

the free Hamiltonian, in the rest frame of matter. The 
scalar field, </>(x), represents quasi-particles whose mass is 
modified by the medium, and propagate in a momentum- 
dependent way. The in-medium mass, m*, is related to 
the vacuum mass, m, via 

m 2 (|k|)=m 2 -<5Af 2 (|k|). 

The mass-shift is assumed to be limited to long wave- 
length collective modes: (5M 2 (|k|) < m 2 if |k| > A s . 
As a consequence, the dispersion relation is modified to 
n 2 , = u! 2 . — <5M 2 (|k|), where f2 k is the frequency of the 
in-medium mode with momentum k. 

The in-medium, thermalized annihilation (creation) 
operator is denoted by 6k (6 k ), whereas the correspond- 
ing asymptotic operator for the observed quantum with 
four-momentum = (£j k ,k), ui 2 . = m 2 + k 2 (wk > 0) is 
denoted by et k (a\j. These operators are related by the 
Bogoliubov transformation, i.e., a kl = c kl 6 kl +sl ki 6^ ki , 
which is equivalent to a squeezing operation. For this 
reason, r k is called mode-dependent squeezing parame- 
ter. The relative and the average pair momentum coor- 
dinates are written as 2 = lui — w 2 , qi, 2 = ki — k2, 
Eij = \{t-Oi and K lj2 = ^(k x +k 2 ). For shorten- 

ing the notation, the squeezed functions are denoted by 



r(i,j,x) = -log[(K? jUlM (x))/(K*»(x)u u (x))] 



u ki (x) + 0J kj (x) 



Q ki (x) + O fej . (x) 



(33) 



is the squeezing parameter. Also, n,j is the den- 
sity distribution, which is taken as the Boltzmann 
limit of the Bose- Einstein distribution, i.e., n-*^(x) ~ 

exp{— [Kl ■' u^(x) — fi(x)]/T(x)}, where the symbol (*) 
implies the use of in- medium mass, whereas it is absent 
if there is no mass-shift. 

In the cases of 7r°7r° or <ficf) correlations, where the bo- 
son is its own anti-particle, the full correlation function 
consists of a HBT part (related to the chaotic amplitude, 
G c (l,2)) together with a BBC portion (related to the 
squeezed amplitude, G s (l, 2)), as shown below 



c 2 (ki,k 2 ; 



= i 



= JV 2 (ki,k 2 ) 

JVi(ki)JVi(k 2 ) 
G c (l,2)| 2 |G S (1,2)| 2 



G C (1,1)G C (2,2) G C (1,1)G C (2,2)- 



(34) 



The invariant single-particle and two-particle momentum 
distributions given by 



G C (M) = G c (fc i) fc i )=JV 1 (k i ) 

G c (l,2) = % /Z3 k 7Ij k 7(a ki ak 2 ), 
G s (l,2) = % /Z3 k 7I3k7(a kl a k2 ). 



w k t (a ki a kl 



(35) 



For an infinite, homogeneous and thermalized medium, 
the part of the full correlation function in Eq. (j34|) . cor- 
responding to the Back-to-Back Correlation is written as 



C, 



bosons 



(k, -k) = l 



|CkS k "k + C-kSl k (n_ k + 1) 

rai(k) ni(-k) 



(36) 

The effects of finite size on BBC are considered af- 
terwards in the text. In Ref.[44], the influence of finite 
emission times is discussed, observing that the BBC in 
this case is suppressed when compared to the instant 
emission. Nevertheless, it was also shown that the max- 
imum of the BBC for each value of |fc|, corresponding to 
G 2 (k, — k), still attained significant magnitude, in spite of 
considering the time suppression. We will see this more 
explicitly later. 



D. Back-to- back correlations for fermions 

Previously, the BBC was shown as a different type 
of correlation between boson-antiboson pairs, occur- 
ring if their masses were shifted. In 2001, T. Csorgo, 
Y. Hama, G. Krein, P. K. Panda and S. S. Padula 
demonstrated that a similar correlation existed between 
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fermion-antifermion pairs [451] . if their masses were mod- 
ified in a thcrmalized medium. In the femtoscopic type 
of correlations, identical bosons have an opposite behav- 
ior as compared to identical fermions, resulting from the 
fact that quantum statistics suppresses the probability 
of observing pairs of identical fermions with nearby mo- 
menta, while it enhances such a probability in the case 
of bosons. However, regarding the Back-to-Back Corre- 
lations resulting from squeezed states, a very different 
situation occurs: fermionic BBC are positive and similar 
in strength to bosonic BBC. Besides, contrary to the the 
femtoscopic correlations, the BBC are unlimited. 

The expressions in the fermion BBC case are similar 
to Eq.® and (|5T|). 



-/Vi(kj) =w ki {ai C4 ak < ) ; iYt(kj) = w k< (a k .a k< ) 
7V 2 (ki,k 2 ) =w kl a; k2 (a ki a k2 ak 2 ak 1 ) • 



(37) 
(38) 



In the above expressions, (O) denotes the expecta- 
tion value of the operator O in the thermalized medium 
and a\a,a\a are, respectively, creation and annihi- 
lation operators of the free baryons and antibaryons 



of mass M and Uk = y M 2 + | fc 2 1 , which are defined 
through the expansion of the baryon field operator as 

= (i/v)EA,v,sKie a A,s + w A',-s L,_£) e<£ - if ; v is 

the volume of the system, u x r and v x , _g are the Dirac 
spinors, where the spin projections are A, A' = 1/2, —1/2. 
The in-medium creation and annihilation operators are 
denoted by b',b,b',b. While the et-quanta are observed 
as asymptotic states, the 6-quanta are the ones thermal- 
ized in the medium. They are related by a fermionic 
Bogoliubov-Valatin transformation, 



5 A',-k 



c k 

ZL s * At 



->A,k 



bl, 



A'.-k 



here c\ = cosri, s\ = sinri, and 

|k 1 |AM(k 1 ) 



tan(2n) 



w(ki) 2 -AfAAf(ki) 



(39) 



(40) 



is the fermionic squeezing parameter. Note that in the 
fermionic case, the squeezing parameter is the coefficient 
of sine and cosine functions, differently than the bosonic 
cases in which appeared their hyperbolic counterparts. 
In Eq. (|39|) A is a 2 x 2 matrix with elements A\ x \ 2 = 
X^cr -kix^ , where ki = ki/|ki|, \ is a Pauli spinor and 
X = —io~ 2 X- Since r is real in the present case, we drop 
the complex-conjugate notation in what follows. 

In order to evaluate the thermal averages above, the 
system is modelled as a globally thermalized gas o 
quasi-particles (quasi-baryons). In this description, the 
medium effects are taken into account through a self- 
energy function, which, for a spin-i particle (we will fo- 
cus on proton and anti-proton pairs), under the influence 
of mean fields in a many-body system, can be written as 
£ = T s +7°S° + 7 4 £\ In this expression, T° is a weakly 
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FIG. 1: The plot shows back-to-back correlations of pp 
(fBBC) and of 0- meson pairs (bBBC), as a function of the 
in-medium modified (p or <f>) mass, m», for \k\ = 800 MeV/c. 
The dependence of the fBBC on the net baryon density is 
shown for three values of ps- In both cases T = 140 MeV 
and At = 2 fm/c. The plots were extracted from Ref. |45j . 



momentum dependent function which, for locally ther- 
malized systems that we are considering, has the role of 
shifting the chemical potential, i.e., /x* = fi — TP. The 
vector part is very small and is neglected. The scalar part 
can be written as T s — AM(k). Within these approxi- 
mations the system can be described with a momentum- 
dependent in-medium mass, m >t (k) = m — AM (|k|). 

We are mainly interested here in the study of the 
squeezed correlation function, which corresponds to con- 
sidering only the joint contribution of the first and third 
terms of the rhs of Eq. (|34|) . In the fermionic case and for 
an infinite, homogeneous thermalized medium, the BBC 
part of the correlation function is written as 



(+-) 



j , j ,.,„ ons (ki,-k 1 ) = l+ [l + (2At«*) 
(1 - n k - n k ) 2 (c k Sk) 2 



c k n k 



(1 



n k )] [c k nk 



+ 4(l-n k )] 



}• 



(41) 



where n k - e x P [ ( n k - M »)/T]+i ! nk ~ ex P [(n k + M »)/T]+i 
in terms of which the net baryonic density is written as 
Pb = (fl/y)Ek { n k^ n k}- In Eq. (JH]) we have included 
a more gradual freeze-out by means of a finite emission 
interval, similarly to what was done in Ref [3] , which has 
the effect of suppressing the BBC signal. 

For a numerical study of the fermionic back-to-back 
correlations, fBBC, we considered, for simplicity, momen- 
tum independent in-medium masses, i.e., m* = M—AM. 
In Fig. (p} we show fBBC for pp pairs as a function of the 
in- medium mass m*, for three values of the net baryonic 
density ps'- for the normal nuclear matter, one tenth of 
this value and for the baryon free region, i.e., ps — 0. We 
show in the same plot results for the bosonic case, bBBC, 
corresponding to cf> meson pair, whose mass is close to the 
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proton mass and was the example used in Ref . [441] . 

We see from Fig.||TJ that fBBC and bBBC are, indeed, 
both positive correlations, with similar shape, and of the 
same order of magnitude. We also observe that fBBC 
is strongly enhanced for decreasing net baryonic density, 
being maximal for ps ~ 0, i.e., for approximately equal 
baryon and anti-baryon densities. 



E. Flow effects on back-to-back correlations 

In the previous discussions, an infinite and homoge- 
neous medium was considered. However, we know that 
the systems produced in high energy collisions, includ- 
ing the ones at RHIC, have finite sizes. Thus, it would 
be important to test if the BBC signal would survive 
when more realistic spatial and dynamical hypotheses 
were considered. For pursuing this purpose, we stud- 
ied the effects on the squeezing parameter and on the 
back-to-back correlation of a finite size medium moving 
with collective velocity (46|. For this, a hydrodynamical 
ensemble was adopted, in which the amplitudes G c and 
G s in Eq. and (|55|) were extended to the special 

form derived by Makhlin and Sinyukov [4p| ] . 



4>4> Back— to — Back Correlation 



G c (l,2) 



1 



(2tt) s 



+ |«_i,_ a | 2 (n-i,_a + l)}, (42) 
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FIG. 2: The effect of a finite emission interval on the back to 
back correlation function, as compared to instant emission, 
is illustrated by the two plots. The dashed curves have been 
reduced by a factor of 400, and the solid curves correspond to 
the suppression by a finite emission duration, of about At ~ 
2 fm/c. The plot in (a) shows this effect in the absence of flow. 
The plot in (b) shows the corresponding result when flow is 
included, with <u>= 0.5; the other parameters adopted to 
produce the curves are R = 7 fm/c, T = 140 MeV. The plots 
were extracted from Ref. 14611 . 



G s (l,2) = -^3 J dV,W^ 2 e M >^{ S « l/2rin _ l!2 
+ ci_2sl 2 ,i(ni-2 + l)}- (43) 

In Eq.(|42]) and d i a^{x) = d 3 W (x; T f ) F{r f )dTf 

is the product of the normal-oriented volume element 
depending parametrically on r/ (the freeze-out hyper- 
surface parameter) and the invariant distribution of that 
parameter F(rf). We consider two possibilities: i) an 
instant freeze-out, corresponding to F(t) — 5(t — To); 
ii) an extended freeze-out, with a finite emission inter- 
val, with F(t) = [6{t - T )/At}e-( T - T0 >/ At . These cases 
lead, after performing the integration in dr in Eq. (f42|) 
and (j43|) with weight (Eij e~ l2Ei ' jT ), respectively to: i) 
(ut+Uj) e-K^+^n,. it) {u i + Uj )[l+[{u i +u j YM 2 ]- 1 l 2 . 

We estimate the geometrical and dynamical effects for 
moderate flow on the BBC in the bosonic case, consider- 
ing the in-medium changes of 0-mesons as illustration. 

For small mass shifts, i.e. ("-"*-> < 1, the flow ef- 
fects on the squeezing parameter are of fourth order, i.e., 
O ( E^^gy \ As a consequence, the flow ef- 

fects on n t j can be neglected, and the factor Cij and Sij 
become flow independent, although they could still de- 
pend on the coordinate r through the shifted mass, m* 
(e.g., through T(x), as in hydrodynamics), which is not 
considered here. 



For the sake of simplicity, and trying to keep the re- 
sults as analytic as possible (for details, see Ref. [4r|). 
we made the hypothesis that the mass-shift was inde- 
pendent on the position within the fireball. We further 
assumed that this last one had a sharp boundary, i.e., 
8m — on the surface, and also the density vanishes 
outside the system volume. The spatial integration in 
Eq. AS)) and (|43p extends over the region where the 
mass-shift is non-vanishing, which is not infinitely large. 
For instance, in relativistic heavy ion collisions is a fi- 
nite region V ~ R 3 ~ (5 — 10) 3 fm 3 . We should keep 
in mind that the vacuum term in the integrand vanishes 
outside the mass-shift region, since it is proportional to 
Si j, which is identically zero in that region. On the other 

hand, the terms proportional to n\j (x) are finite. Being 
so, we can extend the integration in Eq. AS)) and (|4"3"|) 
to infinity and, without much loss of generality, we can 
choose for V a Gaussian profile, exp[— r 2 /(2R 2 )]. This 
study was performed considering two situations: in one 
of them, the mass-shift occurs in the entire system re- 
gion, which is considered as 3-D Gaussian with a circular 
cross-sectional area of radius R = 7 MeV/c at its width. 
The second case considers this volume split in two re- 
gions, the mass-shift occurring only in the internal one, 
with R = 5 MeV/c. Although we adopted several sim- 
plifying hypotheses for reducing the problem complexity 
and treat it analytically, the resulting expressions are still 
intricate, so that a graphical presentation of the results 
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FIG. 3: The plots in this panel are similar to the ones on 
Fig. [2] The main difference is that here we assumed that 
the mass-shift occurred only in a smaller part of the system 
volume. The back-to-back correlation is shown as a function 
of the shifted mass m* on top, and as function of both m, and 
the momentum of each particle (ki = — k2 = k), on bottom. 
The plots in parts (a) and (b) illustrate better the behavior of 
the BBC signal seen in parts (c) and (d), for k — 500 MeV/c 
and for |k| = 1000 MeV/c, respectively. In both cases, the 
dashed curve corresponds to (w) = and the solid curve, to 
(it) = 0.5. In (c), the 3-D plot without flow ((u) = 0) was 
considered, whereas in (d) a radial flow with v = (u)r/R = 0.5 
was included. The plots were extracted from Ref. 46]. 



is more beneficial to the reader. A complete discussion 
with full analytical and numerical results can be found 
in Ref.Hl. 

Therefore, Fig. ^ and Fig. ^ summarize the main 
results of that study. The panel in Fig. ([2]) shows the 
case where the mass shift extends over the entire system 
region. It clearly illustrates the dramatic suppression ef- 
fect of a finite emission duration, since the curves cor- 



responding to instant emission had to be reduced by a 
factor 1/400 in order to fit in the same plots . We can 
also see that the effect of flow, within the approximations 
used, also reduces the signal but in a moderate way. 

For illustrating other interesting features found in this 
preliminary analysis and shown in Fig. ([3]), we chose the 
case with two-regions, as briefly delineated above. It also 
allows for some comparison with the result on Fig. ([2]) . In 
parts (a) and (b) of Fig.©, we can see that, depending 
on the momenta of the back-to-back pairs suffering the 
squeezing correlation, the signal is weaker when flow is 
present (for |k| ~ 1000 MeV/c), almost unaffected by 
flow (for |k| ~ 750 MeV/c), or even slightly stronger in 
the presence of moderate flow (for |k| < 500 MeV/c), 
for the set of parameters chosen in these calculations. If 
flow is absent, however, a monotonically increase can be 
observed with increasing values of the pair momenta, |k|. 
A more complete view of the behavior of the maximum of 
the (jxj) back-to-back correlation can be seen in the parts 
(c) and (d) of Fig.©. It is evident that BBC correlation 
function has a steeper growth with momentum in the 
no flow case, for the same values of the shifted mass, 
to*. The moderate flow picture considered in this study 
still causes the growth of the BBC with |k|, but in a 
considerably smaller rate. With this broader panorama 
in mind, it is easier to understand the plots (a) and (b) 
of Figs. ^ and ([3]). Finally, we see that the strength 
of the signal is directly proportional to the size where 
the mass-shift occurs. Although the values of |k| were 
not the same in those two figures, but already knowing 
that the signal grows with increasing |k|, we see that the 
strength in Fig. (0), where the system size is R = 7 fm, is 
bigger than the corresponding ones in Fig. ([3]), in which 
the squeezing region has R s = 5 fm, even for |k| = 1000 
GeV/c. 

Naturally, the above study was mainly a first step to- 
wards better understanding the nature and conditions 
of survival of the squeezing correlation, since it involved 
several approximations in its derivation. In particular, 
it was supposed that the squeezing occurred homoge- 
neously throughout the system region, and independently 
on the particle momentum. Next steps will require a 
modelling instead of this assumption, which will provide 
the volume dependence of the squeezing in a more realis- 
tic way. Besides, we have presented merely the behavior 
of the maximum of the Back-to-Back Correlation func- 
tion. Is is under investigation how the signal behaves for 
finite systems and small values of the average pair mo- 
menta, which hopefully will furnish the optimized form of 
the BBC signal that should be looked for experimentally. 



F. Chiral dynamics and back-to-back correlations 

This sub-section is based on ref. [47} • The environ- 
ment generated in the mid-rapidity region of a high- 
energy nuclear collision might endow the pionic degrees 
of freedom with a time-dependent effective mass. This 
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implies two-mode squeezing, but with a time-dependent 
squeezing parameter. Its specific evolution provides a 
mechanism for the production of back-to-back charge- 
conjugate pairs of soft pions which may present an ob- 
servable signal of the non-equilibrium dynamics of the 
chiral order parameter. The suddenness of the transition 
to the asymptotic quanta is a condition that is released 
in this approach. The important point of the numeri- 
cal investigations along such a model assumptions were 
that suddenness is not a mandatory requirement. The 
BBC signal seems to be strong enough to survive even 
the time-dependence of the effective mass, caused by chi- 
ral dynamics, resulting in an adiabatic mass variation, 
modelled by mass oscillations within an exponentially de- 
creasing envelope. 



G. Experimental constrains on chiral dynamics as 
a suggested explanation of the RHIC HBT puzzle 

Recently, Cramer, Miller and collaborators published 
an interesting work, that included a relativistic quantum 
mechanical treatement of opacity and refractive effects, 
that allowed for a reproduction of STAR two-pion (HBT) 
correlation data and pion spectra in Au+Au collisions at 
y/ snn — 200 GeV colliding energies at RHIC. This in- 
vestigation suggested, that an attractive, real part of the 
optical potential is the critical element needed to repro- 
duce the transverse mass dependence of the sidewards 
and the outwards HBT radii at RHIC. This optical po- 
tential represents the strength of the interactions between 
the pions and the medium. Chiral dynamics in the model 
leads to a temperature and density dependent in-medium 
modification of the pion (pole as well as screening) mass. 
The results were found to be consistent with a system, 
that had a restored chiral symmetry [igj . 

Cramer, Miller and collaborators have suggested to 
check experimentally various new phenomena to see if 
their explanation is correct: they have argued that a pi- 
onic version of the Ramsauer resonances lead to peaks in 
the HBT radius parameters R and R s in the low mo- 
mentum region, 15 MeV < pt < 65 MeV, and have also 
predicted a peaking behaviour in the same region in the 
transverse momentum spectrum of pions [H, [3 H3| • 

Apparently, there are various other possible tests that 
can be performed to see if this connection of pion HBT 
radii to chiral dynamics is a unique, and fully correct 
explanation, or not. For example, if in medium mass 
modifications of pions is the reason for the effect, which 
depends on the interactions of pions with the medium, 
then similar interactions but of different strength should 
be present between kaons and the medium, hence kaon 
and pion HBT radii are not expected to show a similar 
transverse mass dependence. On the other hand, if the 
HBT radii are equal due to asymptotic properties of ex- 
act solutions of relativistic hydrodynamics, as proposed 
by the Buda-Lund hydro model [17| . kaons, pions, and 
all other heavier particles should show the same effective, 



and transverse mass dependent, approximately spheri- 
cally symmetric source size at large transverse masses, 
{Rout — Rside — Riong oc l/^/mt). So the calculation of 
the transverse momentum dependence of the pion HBT 
radii above the presently limited to 600 MeV upper limit 
in transverse momentum from this model, and the com- 
parison of this value to the PHENIX results on these radii 
up to 1.2 GeV in transverse momentum, and a theoretical 
and experimental combined effort to do a similar study 
for kaons could help to clarify the case. 

However, such an investigation inevitably needs a co- 
operation between theorists and experimentalists, and 
any possible conclusions will have a definite model de- 
pendence. From the experimental point of view, it is 
better to find a clear-cut effect, which has to exist, if in- 
deed chiral dynamics is the key element in the success of 
the Cramer-Miller approach. As should be clear from our 
presented review on squeezed correlations, back-to-back 
correlations of particle- antiparticle pairs were shown to 
exist in all scenarios that discuss in-medium mass modi- 
fication of hadrons, regardless of the details of the calcu- 
lation (bosons or fermions, expansion, finite duration of 
particle emission, in-medium mass modification only in a 
fraction of the total volume, time dependent in-medium 
mass modification due to chiral dynamics). Thus the ex- 
perimental determination of the in-medium mass modifi- 
cation with the help of back-to-back particle anti-particle 
correlations could give a clear-cut answer, if the success of 
the Cramer-Miller model is indeed due to a chiral phase 
transition, or perhaps due to other, less unique elements 
of their model. 



VII. SUMMARY: A NEW CORRELATION 
SIGNATURE OF A FREEZE-OUT FROM A 
SUPERCOOLED QUARK-GLUON PLASMA 

Currently available correlation data seem to exclude 
the possibility of a strong first order phase transition in 
Au+Au collisions at RHIC and in Pb+Pb collisions at 
top CERN SPS energies. The proposed signature of a 
possible second order QCD phase transition has not yet 
been investigated in detail, the excitation function of the 
Levy index of stability needs to be determined experi- 
mentally. From lattice QCD calculations there is an in- 
dication that, in case such a critical point is reached at a 
certain value of the center of mass energy then, above this 
point, a cross-over type of transition should be expected. 
Presently, a cross-over type of equilibrium hadroniza- 
tion, or, a sudden non-equilibrium hadronization mecha- 
nism are both possible scenarios for hadron production in 
ultra-relativistic heavy ion collisions at RHIC. We have 
focused here on the latter possibility, and discussed a 
quantum optical characteristics of such a scenario, argu- 
ing that it leads to vanishing in-medium mass modifica- 
tions of hadrons. As a consequence, the back-to-back cor- 
related particle- antiparticle pairs should disappear at the 
onset of such a sudden deconfincmcnt transition, charac- 
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teristic of quark-gloun plasma hadronization and freeze- 
out from a deeply supercooled, negative pressure state. 
These back-to-back correlated particle-antiparticle pairs 
should on the other hand appear, if a chiral phase tran- 
sition is present in y/sJfN — 200 GeV Au +Au collisions 
at RHIC. Thus they can also be utilized as model in- 
dependent experimental tools to test the validity of the 
Cramer-Miller approach. 
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